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Abstract 



We study the influence of Gribov copies in the Coulomb gauge on 
the smeared hadronic correlation functions that are involved in the de- 
termination of the B meson decay constant. We find that the residual 
gauge freedom associated to Gribov copies induces observable noise ef- 
fects, though at the level of numerical accuracy of our simulation these 
effects are not relevant to the final determination of fs- Our results 
indicate that such effects may become important on bigger lattices. 
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In the last few years numerical studies of lattice gauge-fixing ambiguities 
have been performed by several groups [1-4]. Such studies are interesting 
for many reasons. On the theoretical side, the existence of these ambiguities 
parallels an analogous problem in the continuum formulation of nonabelian 
gauge theories [5-10], and the fact that these nonperturbative features of the 
theory are reproduced both by continuum and lattice regularized models is 
in itself reassuring. Moreover, in the analytical study of the continuum limit 



of lattice gauge theories ||ri|, |T2[ gauge-fixing is inescapable and the study of 
Gribov ambiguities turns out to be an important point in such a program. 
Finally, we mention that gauge dependent quark and gluon matrix elements 
may be used to derive renormalization conditions, as recently pointed out in 
[|T^ , and gauge- fixing ambiguities may play a role in this kind of application. 

From a more practical point of view, we observe that our current under- 
standing of QCD phenomenology on the lattice very often proceeds through 
Monte Carlo simulations involving gauge dependent operators, like hadron 
wavefunctions and, in particular, "smeared" sources for hadronic correla- 



tion functions |T^, |I5|. In such cases, if the gauge condition we implement 
numerically (e.g. the Coulomb gauge) does not correspond to a complete 
gauge-fixing, in the sense that the gauge-fixing algorithm may converge ran- 
domly to any configuration in a set of Gribov copies, then the value of the 
operators will depend on which copy gets selected by the algorithm. As we 
will discuss below, such residual gauge freedom acts as a source of statisti- 
cal noise in the Monte Carlo average of those physical quantities that are 
extracted from gauge dependent quantities. 

The physics of the B meson on the lattice is a suitable sector for an 
investigation of the above issue. The lattice calculation of the decay constant 
fs is performed in the static approximation, in which the heavy constituent 
quark has infinite mass and, consequently, the heavy quark propagator is 



given by a product of time- like links [16] 



Given the local, gauge invariant "heavy-light" axial current 



Af^{x,t) = b{x,t)-ff,-f5q{x,t) (1) 

(where b, q indicate field operators associated to a heavy and a light quark 
respectively) one is interested in the asymptotic behaviour in time of the 
2-point correlation function at zero spatial momentum: 

C^^{t) = Y.< 0\T{A^ix, t) A^{0, 1)}|0 >^ Z^Z^ e-^^* (2) 
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In principle Jb can be obtained from Z^, but in actual simulations it has not 
been possible to isolate the lightest pseudoscalar state in the above correlation 
function, thus no direct evaluation of from C^^ can be obtained [IT? 



This problem has been overcome by evaluating on the lattice the correlation 
function associated with the "smeared", gauge dependent, axial current 

A'^ix, t) = ^J2 t)^,l5qix, t)} (3) 

i 

The sum runs over lattice points in a 3-dimensional spatial cube cen- 
tred on {x,t), n being an odd integer. The 2-point function C^^{t) of the 
extended current (^, defined analogously to (0), and the mixed 2-point func- 
tion C^^{t) are characterized by a much better numerical signal. Therefore 
their (asymptotic) large t behaviour 

C^^it) ^ Z^Z^ e"^^* and C^^{t) ^ Z^ Z^ e'^^' (4) 

makes it possible, by taking ratios of these gauge dependent correlation func- 
tions, to extract the gauge invariant Z'^ITSl . 



Such calculations are typically performed in the lattice Coulomb gauge. 
In fact, it is well known that if one tries the numerical evaluation of C^^ and 
C^^ without previously fixing the gauge, one experiences large fluctuations 
induced by the gauge freedom, with no detectable signal. It is therefore con- 
ceivable that the residual gauge freedom associated to an incomplete gauge 
fixing prescription (as we will see the lattice Coulomb gauge, containing a 
large number of Gribov copies, is a good example) would still induce un- 
wanted fluctuations in the gauge dependent correlation functions C^^ and 
C^^ and in their ratios, thus contributing to the overall error bars for Z^ 
and fs- 

The aim of this letter is to test numerically this possibility by evaluating 
C^^ and C^^ in the Coulomb gauge but on different Gribov copies of the same 
thermalised configuration. We will follow closely the procedures adopted in 



ref. [^. Before focusing on this, it is necessary to discuss the method that 
we have adopted for the generation of such copies. 

We consider a lattice with spatial volume Vg and T sites in the time 



direction. The standard lattice implementation of the Coulomb gauge [19 



|l[] consists of taking a thermalised link configuration {^7} and iteratively 
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applying gauge transformations to it in order to minimize the quantity 



and the gauge field transforms as Uf{x,t) = g{x,t)Ui{x,t)g\x + i,t). The 
lattice Coulomb gauge condition is automatically satisfied when the trans- 
formed lattice {U^} is such that F[U^], considered as a function of the gauge 
transformations g, attains a local minimum. 

We emphasize that by deciding to implement the Coulomb gauge through 
the minimization of F[U^] one has already eliminated some of the lattice 
Gribov copies of such a gauge, since any stationary point of F[U^] (including 
saddle points and local maxima) would be a configuration satisfying the 
gauge condition, hence a Gribov copy. For this reason the minimization of 
F implements a gauge containing additional constraints with respect to the 
standard Coulomb gauge. This corresponds to the original proposal made 
by Gribov for the quantization of the theory in the continuum, but it still 
turns out in both lattice and continuum models that one has not eliminated 
all the gauge freedom, since in general F[U^] as a function of g has many 
local minimafl, These are the Gribov copies we will be dealing with. 

A set of Gribov copies for each Monte Carlo thermalised configuration 
can be generated in several ways. One may perform random gauge transfor- 
mations on the thermalised configuration before fixing the gauge, or change 
the route of updating the lattice sites when gauge fixing, or even vary a 
relevant parameter of the gauge fixing algorithm, as in the overrelaxation 
procedure |pO|, 

Eqs. and (|^) imply that, given a starting link configuration, the 
Coulomb gauge condition is implemented independently on each timeslice. 
This is because the links in the time direction which connect adjacent times- 
lices do not appear in the definitions of F[U^] and f[U^]{t). In other words, 
each timeslice of a given configuration is endowed with its own scenario of 
Gribov copies, so that the pattern of their occurrence in the Coulomb gauge 
is richer than in the Landau gauge, in which they are defined globally on the 
whole lattice. 



t=l 



(5) 



where 
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Consider then a gauge transformation g, that transforms the original 
Monte Carlo configuration {U} into the Coulomb gauge- fixed one {U^}. 
Since Coulomb gauge-fixing is an independent process on each timeslice, it is 
convenient to think of the set of matrices ^ as a union of subsets ^ = IJ^ gt, 
where the gt implement the gauge on each timeslice 0. It is clear that if 
9^ = Ut9h 9^ = Ut gh 9^ = Ut 9t are distinct gauge trans- 

formations, that rotate the configuration {[/} into different realizations of 
the Coulomb gauge, then any combination g'^°™-^ = Ut 9ti where for each 
value of t the index i can take any value between 1 and A^, is a lattice gauge 
transformation that implements a new distinct Coulomb gauge. In other 
words, a Coulomb gauge transformation for the entire lattice can be built 
up from any of the combinations of partial, fixed-time Coulomb gauges. 
This property allows us to single out of the possible choices two "spe- 
cial" realizations of the Coulomb gauge as follows: for each timeslice t of the 
configuration {[/} we select from the set of fixed time Coulomb gauges the 
transformation ^J"*" such that f\U^™"]{t) takes the smallest value. In such 
a way one can define g^^^ = Ut ^r*" Coulomb gauge transformation 

such that takes the smallest possible value. Analogously, one can 

build up from our set of fixed time gauge transformations defined as 

the gauge transformation such that takes the largest possible value. 

In the following we will analyze the variation of the gauge dependent 
operator C^^(t) when evaluated on {[79""""} and {U^""'"}, comparing it to 
the corresponding variation of f{t). In this way we can estimate the order of 
magnitude of the fluctuations of C^^(t) induced in an actual simulation by 
the residual gauge freedom associated to the Gribov copies. Such an estimate 
is based on the assumption that the configurations {[/a""'"'} and {U^""'"}, that 
have been constructed so that they maximize A[F] = — 
also maximize the variation of C^^{t). In some sense, we are assuming that 
F is a reliable measure of the "distance" between Gribov copies. In any case, 
our estimate should provide a lower bound for the magnitude of the effect. 
We have also obtained numerical results for C^^, that are qualitatively of 
the same type as those for C^^ but characterized by a worse signal, so that 
we prefer to focus the discussion on C^^. 

We have considered five SU (3) lattice configurations generated on a lat- 
tice of size 10^ X 20 at /3 = 6.0. The lattice has been doubled in the x and t 
directions before evaluating the quenched light quark propagator for Wilson 
fermions ai K = 0.1515. For these lattice parameters, the best smearing 
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size is 72 = 7, according to |T8|]. We have fixed the gauge to an accuracy 
< 10~^° for each time slice; this eliminates the possibihty of misin- 
terpreting fluctuations of / related to the poor level of gauge fixing as the 
evidence of distinct Gribov copies. 

In the Tables (|I|) and (|^) we report numerical results (half lattice) from 
two of the five configurations that we have analyzed. As can be seen from the 
second column of the tables, up to 6 distinct values of f{t), i.e. 6 different 
Gribov copies, have been obtained on some timeslices. On other timeslices it 
has not been possible to obtain as many as 6 distinct values, not even after 
generating 300 random configurations, and in some cases no configurations 
with different values of f{t) were found. This last feature takes place in all 
5 thermalised configurations, but obviously not on the same timeslices. 

The two tables show different behaviors of C^^{t): in Table (|I]) we see 
that C^''^{t) fluctuates only when the value of f{t) fluctuates, while in Table 
(0) C^^{t) always changes, even at those times when f{t) does not. The 
difference is due to the fact that, contrary to the case of Table (|I|), in Table 
(0) / fluctuates on the flrst timeslice, so that C^^{t) is expected to change, 
since one of the smeared sources is always located on the flrst timeshce. 

In particular, the results of Table confirm the existence of a correspon- 
dence between the fluctuations of f{t) and those of the smeared hadronic 
sources at the same t, that generate the fluctuations that we measure on 
C^^{t). This justifles a posteriori our decision to probe the fluctuations of 
hadronic correlation functions by evaluating them on those conflgurations 
that maximize the fluctuations of f{t). 

In all cases examined we found, as expected, that C^^ fluctuates per- 
centually much more than /. In particular, by inspection of Table 1 it turns 
out that while ^^jp is always less than 1%, the percentual variation of C^^{t) 
increases with t, i.e. with the separation of the smeared sources, starting 
from less than 1% for t = 1, going to more than 10% at t = 9 and being 
~ 50% or more for t > 15. 

This is related to the fact that f{t) is the average on the timeslice t 
of a local function of the link variables, hence is sensitive to the "noise" 
effect induced by the Gribov copies on that single timeslice; on the contrary 

CSS 

is a 2-point correlation function of fermionic operators and picks up 
a "noise" effect from the Gribov copies on both the timeslices where the 
smeared sources are located. Moreover, while the signal that we want to 
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extract from C^^{t) dies out exponentially with t, the "gauge noise" effect 
is expected not to depend on the separation of the sources, thus becoming 
more and more important percentually at large separations. 

However, the gauge noise does not produce on our present lattice a sig- 
nificant difference in the physics that we are considering, since the signal is 
very noisy at large times, thus only the time interval t = 5 — 10 may be 
used for the determination of and Z^, as in ref.[|T8l. In such interval the 



gauge noise effect is not the dominant contribution to the numerical uncer- 
tainty contained in the procedure for the determination of fs, that is mostly 
affected by errors due to the fitting of the correlation functions. 

Nonetheless, it is worth noting that the measurement of fs is not spoiled 
by the gauge noise also because in the range of interest for t, that is between 5 
and 10, in many cases we have not been able to find copies. This fact reduces 
drastically our capability to fully evaluate the gauge noise effect expected on 
this lattice. 

On the basis of our results, we conclude that the residual gauge freedom 
in the standard lattice implementation of the Coulomb gauge indeed shows 
up as a noise effect on smeared hadronic correlation functions. 

While the gauge noise is not too important in the framework of the ac- 
tual measurement that we have discussed, since this is plagued by many 
sources of error, on the other hand, for the reasons that we have sketched 
above, such an effect might sensibly affect measurements on bigger lattices. 
In fact, increasing the lattice size would typically reduce the standard sta- 
tistical noise, allowing in principle to perform measurements from large time 
separations, but then the gauge noise effect (that is not expected to disap- 
pear when approaching the continuum limit) may become a major source 
of fluctuations and provide a relevant contribution to the error bars on the 
physical quantities evaluated from smeared correlations. We plan to further 
investigate the subject by increasing the lattice size and the statistics and 
also by considering other smeared operators. 

We thank V. Lubicz, G. Martinelli and M. Testa for helpful discussions. 
We also thank G. Martinelli for allowing us to use an analysis program. 
C.P. acknowledges financial support from C.N.R. and I. N.F.N. 
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Table 1: The values per timeslice of the gauge dependent quantities f{f) 
and C^^{t) for the two Gribov copies f/^""" and ]\[ jg number of 

distinct values for f{t) we found on each timeslice. 
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Table 2: Same as in Table 1 for a different thermalised configuration. 
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